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ABSTRACT 



We derive a bound on the total number of negative energy bound states in a potential in 
two spatial dimensions by using an adaptation of the Schwinger method to derive the Birman- 
Schwinger bound in three dimensions. Specifically, counting the number of bound states in a 
potential gV for g = 1 is replaced by counting the number of g^s for which zero energy bound 
states exist, and then the kernel of the integral equation for the zero-energy wave functon is 
symmetrized. One of the keys of the solution is the replacement of an inhomogeneous integral 
equation by a homogeneous integral equation. 
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1 Introduction 



In a previous paper [1], K. Chadan, N.N. Khuri and ourselves (A.M. and T.T.W.) obtained a 
bound on the number of bound states in a two-dimensional central potential. This bound has 
the merit that, for a potential gV, the coupling constant dependence for large g is optimal, i.e., 
the same as the one of the semiclassical estimate [2]. Previous work on the subject was done 
by Newton [3] and Seto [4]. We also obtained a bound for the non-central case, but only by 
using a rather brutal method which consists of replacing the potential by a central potential 
which is defined, after choosing a certain origin, by 

V c (r) = InfV(x) 

\x\ = r (1) 

Because of the monotonicity of the energy levels with respect to the potential, putting V c in 
our formulae will give a bound for the potential V. However, for potentials with singularities 
outside the given origin, this may lead to no bound at all. Our attention has been attracted 
by the fact that in condensed matter physics problems exist, where counting the bound states 
on a surface may be useful [5], but where it is very unlikely that the potential will be central, 
even approximately. 

In the present paper, we obtain a bound on the number of bound states in a non-central 
two-dimensional potential, using an adaptation of the Schwinger method to derive the Birman- 
Schwinger bound [6] in the three-dimensional case. The condition under which we obtain a 
bound is 

J jd 2 xd 2 y\V(x)\(en\x-y\) 2 \V(y)\<oo. (2) 

This condition is non-linear , just like that of Birman and Schwinger, but we show in the 
Appendix that it follows from the linear conditions 

J d 2 x (£n(2 + \x\)) 2 \V(x)\ < oo 
and 

J d 2 x V R {\x\)£n-\x\ < oo, (3) 

where Vr is the circular decreasing rearrangement of |V| (for the definition of Vr, see the 
Appendix). 

Condition (2) has already been proposed by Sabatier [7]. Condition (3) appears in a forth- 
coming work by N.N. Khuri, A. Martin, P. Sabatier and T.T. Wu, dealing with the scattering 
problem. It has the advantage of showing more clearly what kind of behaviour the potential is 
allowed to have at short and large distances. 

The strategy of Schwinger consists of counting the number of zero-energy bound states for 
a potential gV for < g < 1 instead of the actual number of negative energy bound states 
for g — 1. In three dimensions these two numbers are equal. Indeed, let < gi, g2, ■ ■ ■ g n < 1 
be the coupling constants for which we have zero-energy bound states. Each g^ is the origin 
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of a bound state trajectory in the E — g plane, E^g), with E^gj) = 0. These trajectories are 
monotonous decreasing: 

f=/W, (4) 
by the Feynman-Hellmann theorem, but 

g t J Vi) 2 d n x = E - J \Vi/j\ 2 d n x < for E < 0. (5) 

This shows that the number of negative-energy bound states is exactly the same as the 
number of g^s < 1. At the crossing of any pair of trajectories there is no problem because of 
their monotonicity. 

The same result holds in two dimensions with one modification: any attractive potential 
(i.e., / d 2 xV < 0), has a bound state for arbitrarily small g, with a binding energy going to 
zero for g — > like — exp — (C/g) [1]. At E — it disappears and is not included in Schwinger's 
accounting, so we have to add one unit. 

Since we only want a bound on the number of bound states, we can always replace V(y) by 
-V-(y): 

V~(y) = 0for\/>0 

V~(y) = -V(y)iorV<0 (6) 
Using — |V(y)| instead of —V~{y) gives a more crude bound. 

It can be shown that the general solution of the zero-energy Schrodinger equation 

-Aip - V~ip = (7) 

in the equivalent integral form 

i>{x) =C-^- J £nk \x - y\V-(y)^(y)d 2 y (8) 
has a general asymptotic behaviour for \x\ — > oo 

iP(x) - C ~ -£nk \x\ i- J V-(y)iP(y)d 2 y + o(l) (9) 

under the condition 

J V~(y) (£n(2 + \y\)) 2 d 2 y < oo (10) 

which follows from condition (1) as shown in the Appendix. 

Zero-energy bound states are characterized by the fact that ip is bounded . Hence we get 
the necessary condition: 

J V-{y)^{y)d 2 y = Q (11) 

Now we have two possibilities: 
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• I. At infinity if){x) — > and hence C = and those bound states wave functions satisfy a 
homogeneous integral equation 

*l>i(x) = ~^J d "y in (ko\x - y\) V-(y)Mv) (12) 

(notice that the scale factor k disappears because of condition (11)). 

This is what happens in the case of a central potential for a non-zero azimuthal angular 
momentum m. 

• II. At infinity, ip(x) — > C, with C ^ 0. In this case, the bound state wave functions satisfy 
an inhomogeneous integral equation. This case has been described in Ref. [8], where it 
is shown that for a central potential in two dimensions, the m = phase shift has the 
universal behaviour 

except if there is a zero-energy bound state of type II. Then 

5(k)£nk -> 0. (14) 

In Ref. [8], a much stronger result is stated. This much stronger result, however, holds 
only for a very rapidly decreasing potential. 



2 Counting Bound States in Case I 

Following Schwinger, we symmetrize the kernel of the integral equation: 

<t>i{x) = 9i J K(x,y)(f)i(y)d 2 y (15) 

with 

(f)i(x) = \/V~(x) ipi(x) 
K(x,y) = -i- y/v~(x) £nk \x-y\ y/V-(y). (16) 

If V~(x) vanishes in some regions, it seems impossible to go back from <pi to ipi- However, this 
can be remedied by defining 

V~(x) = V~(x) +eexp-//|x|. (17) 

Since the bounds we shall get are continuous in V, we can take the limit e — > at the end. 
K can be written as 

tf = E-|&X&|+i2 (18) 

9i 
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R is a sum over states which do not satisfy (11). The 0j's in themselves do not form a complete 
set. If we define a by 

a(x) = , V (19) 
y/S V-(y)<Py 

we have 

< a\4)i >= (20) 
from property (11), and naturally < a\a >= 1. 

If we define Tr, a trace restricted to the 0j's, we have 

TrK = £ - > ^ — > iVj, 

TV/ being the number of bound states of type I. However, this trace turns out to be divergent 
because of the logarithmic singularity of the kernel (the same happened in Schwinger's original 
work!), and we follow Schwinger to iterate the integral equation (15): 



M x ) = 9i j K(x, y)K(y, z) <j)i(z)d 2 yd 



2 z 



and then 
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TrK 2 = E > JVj (21) 
9i 

Forgetting the "hat" on the trace still gives a bound because K 2 is a positive operator (contrary 
to K!), but this bound depends on the scale parameter fc entering in the logarithm. Among the 
missing states in Tr is the state \a >, orthogonal to the 0j's, and this one should be removed 
from the complete trace. In this way, we get 

Nj < TrK 2 - < a\K 2 \a >, 

or more explicitly 
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^ J V-(x)(£nk \x-y\) 2 V-(y)d 2 xd 2 y 
(27r)2 j V -l x)d 2 x j V-(x)(lnk \x-z\)V-(z)(lnk \z-y\)V-(y)d 2 xd 2 yd 2 z(22) 



It is visible that the second term is negative as we announced. 
Rewriting Nj as 



. . / d 2 xd 2 yd 2 zV (x)V (y)V (z) (in k \x — y\) 2 — in k$\x — z\£n k \y — z\ 
[z)d z z J L 



JV-(z)<Pz 

we see that (22) is manifestly independent of the scale factor k . 
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3 Counting Bound States in Case II 



At first it would seem that Schwinger's technique will not work because, in Eq. (8), the constant 
is not zero and therefore we deal with an inhomogeneous integral equation which can be written, 
after the same changes of variables as in Section 2, given by (16) and (19): 

(pi = Qa + giK(pi, (23) 

with, again, 

< a\4>i >= and < a\a >= 1 (24) 

Equation (24) is precisely the key property which will make it possible to replace (23) by a 
homogeneous equation. 

Again, the 0;'s corresponding to different g^s are orthogonal because 

< (fiilcfij >= Cj < <pi\a > +gj < <t>i\C<t>j > 
= Cj < 4>j\a > +gi < 0i|C0j > . 

Hence, from (24): 

(---]< <f>i\<l>3 >= 0- (25) 
\9i 9jJ 

Let us call S the Hilbert space associated to the integral equation (23), and construct a new 
Hilbert space by removing the element a: 

S = S'(B {a} (26) 

We want to define a new operator K' acting in S'. Let 

b = Ka. (27) 

Notice that 

C C 
<b\(f) i >= t -< a \a>= -. (28) 

9i ' 9i 

We try 

K' — K — \b >< a\ — \a >< b\ + C\a >< a\ 
where C will be chosen so that 

K'a = 0. (29) 

We have 

K'\a >= \b> -\b> - < b\a > \a > +C\a >, 

and hence we take 

C =< b\a >=< a\K\a > (30) 
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K' is Hermitian like K, and we get 



giK'lfa > = gi\K\<f>i > -g t < b\<f>i > \a > 
= gi\K\<pi > +d\a > . 

Hence 

giK'lfa >= |0, > (31) 

which is homogeneous. 

To get a bound on the number of bound states of type II, we have to get a bound on trace 
K' 2 (not surprisingly, trace K' is divergent). It is a lengthy but straightforward exercise to 
calculate that trace, which gives 

N n < tiK 2 - 2 < a\K 2 \a > +(< a\K\a >f (32) 

The last two terms give an overall negative contribution. The first term is the same as the one 
appearing in Nj. It is easy to see that the right-hand side of (32) is independent of the scale 
parameter k entering into the kernel K. Finally, let us notice that the treatment of case II 
contains case I because, in the argument, it has never been said that Cj ^ 0. Equation (31) 
holds irrespective of whether Cj = or C; ^ 0. Notice that the bound on Nj is larger than the 
bound on N n . Therefore the bound on Nj becomes completely obsolete. 



4 Concluding Remarks 

If we include the bound state with evanescent energy for zero coupling constant, we get the 
bound 

N <l + trK 2 - 2 < a\K 2 \a > +(< a\K\a >) 2 . 

Dropping the last two terms still gives a scale-dependent bound - which can be minimized with 
respect to the scale - which precisely appears in condition (2), itself following from the linear 
condition (3) as shown in the Appendix. 

Conditions (2) and (3) both allow a potential behaving like 

1 

r 2 (£nr)^ 

at infinity, with local singularities not worse than 

1 

\ r — r Q \ 2 (\£n\r — r | |) 2+e ' 

e positive, arbitrarily small. Both conditions are violated for e < 0. However, we shall see in 
the Appendix that (2) is definitely weaker than (3). 

Our bound has the merit of being valid for the non-central case, which, as we said in the 
Introduction, is important for solid-state physics. However, for a potential gV, it behaves like 
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g 2 for large g, while in Ref. [1], in the central case, we get a bound behaving like g. In Ref. [1] 
we make a conjecture which is very far from being proved, but clever mathematical physicists 
might prove it or something similar. The present work should be considered only as a first step 
which could possibly give reasonable results for not too large g. 
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Appendix: Comparison of condition (2) and condition (3) 

Condition (2) is 

/ = J d 2 xd 2 yV-(x)(en\x-y\) 2 V~(y) < oo (A. 1) 

Condition (3) is a set of two conditions: 

J d 2 x (£n(2 + |:r|)) 2 V~{x) < oo (A. 2) 

J d 2 xV R {\x\)£n-{\x\) < oo. (A. 3) 

In (A. 3) we use: 

-£n-(\x\) = Oif|x|>l 

= — £n\x\ if \x\ < 1. 

— Vft(|x|), the circular decreasing rearrangement of V~(x). Since this notion is not very well 
known among physicists, let us remind the reader that Vr(|x|) is a decreasing function of 
such that 

fjL(V R (\x\)>A) = n(V-(x)>A),WA, 

where fi is the Lebesque measure. In more familiar terms, the rearranged Mont Blanc would 
be a mountain with axial symmetry, with a single peak, such that the surface between the level 
lines would be the same as the surface between the level lines of the original Mont Blanc (rather 
awfully dull!). 

We shall prove first that the convergence of / in (A.l) follows from the convergence of (A. 2) 
and (A. 3). More exactly, we shall get an explicit bound on (A.l) in terms of (A. 2) and (A. 3). 
We write 

/ = /+ + /_ (A. 4) 
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with 



/+ = Jd 2 xd 2 yV~(x) (£n + \x -y\) 2 V~(y) 
/_ = J d 2 xd 2 yV-(x) (£n~\x -y\) 2 V-(y) 



(A. 5) 
(A. 6) 



£n has already been defined. £n + (a) = £na for a > 1, = for a < 1. It is elementary to get a 
bound on I + from (A. 2) only. Indeed, 



< £n + \x -y\< £n + (\x\ + |y|) < £n{2 + \x\) + £n(2 + \y\) 

(ln + \x - y{f < 2 [(£n2 + \x\) 2 + (£n(2 + \y\)f 

I + <Aj d 2 xV~(x) J d 2 yV-(y) (£n(2 + \y\)) 2 . 
The convergence of the right-hand side of (A. 8) follows directly from (A. 2). 



and thus 
Hence 



(A. 7) 
(A. 8) 



Concerning J_, we use a rearrangement inequality due to Luttinger and Friedberg [9], which 
says 

J J A{x) B(\x - y\) C(y)d 2 xd 2 y < J A R {\x\)B R {\x - y\) C R (\y\) d 2 xd 2 y (A. 9) 

where A, B, C are non-negative functions and A R , B R , C R are their decreasing rearrangements. 
Since £n~ and (£n~) 2 are decreasing functions of their argument they are their own rearrange- 
ment. Hence 

/_ < J d 2 xd 2 y V R (\x\) (£n'(\x - y\)f V R (\y\) (A. 10) 

In (A. 10), we can carry out first the angular integration, the angle (x, y) appearing only 
in £n~ . However, to be able to do that easily we have to sacrifice some information, i.e., use 
(£n~\x — y\) 2 < (£n(\x — y\)) 2 - We have to calculate 



r dB_ 
J 2^ 



x 



\y\ 2 — 2\x\ \y\ cos 9 



We have 



£n(\x\ 2 + \y\ 2 - 2\x\\y\ cos9) = £n(\x\ - \y\e ie ) + £n(\x\ - \y\e~ 



Assume \x\ > \y\. Then we get 



£n{\x\ 2 + \y\ 2 - 2\x\\y\ cos9) = 2 



£n\x\ — E 



\x\ 



cos n9 



n 



(A. 11) 



Hence, if \x\ > \y\, using the orthogonality of the cos n9: 



/ — (£n(\x\ 2 + \y\ 2 - 2\x\\y\ cos9)) 2 = A(£n\x\) 2 + 2 V f M I A 

J 27T n=l\\ X \J U ' 
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We see a dilogarithm, or Spence function, appearing on the right-hand side. However, we only 
need to notice that 

/d6 2 00 1 7T 2 

— (€„|x| 2 + |y| 2 -2|y||x|oosfl) = 4(£n\x\) 2 + 2 ]T - = 4(in\x\) 2 + —. (A. 12) 

n=l n 13 

In this way we get 



/_<(2tt) 2 x2 / NdlxlHdM^dxD^dyl) 

^|x|>|j/| 

Again, we split the integral into 



7T 



4(en\x\y + — 



32tt 2 



/ |x|d|x||y|d|y|^d*l)^M)(^>l) 2 

•/|x|>|j,| 



+ 32tt 2 / \x\d\x\\y\d\y\V R (\x\)V R (\y\)(in+\x\) 2 

J\x\>\y\ 

2 



1 47rl 
^ 3 



I IxIdlxlFfldxl) 



(A. 13) 



In the first term of (A. 13) we can replace (in (\x\)) 2 by in \x\ in \y\, since \x\ > \y\ and since 
in~ is decreasing. 

In the second term, we can drop the restriction \x\ > |y| and notice that 

J d 2 xV R (\x\) (in + (\x\)f < J d 2 xV-(x) (in + (\x\)f . 

Indeed, / d 2 xA R (x)4>(\x\), where </>(\x\) is increasing , is less than / d 2 x A(x) 4>(\x\). Suppose 
that 0(|x|) — > L. Then 

J d 2 xA R {x)(t>{\x\) = J d 2 xA R (x)L- J d 2 xA R (x)(L-(f>(\x\)). 

L — <f>(\x\) is its own decreasing rearrangement and following the well-known properties 

J A R B R d 2 x > J A(x)B(x)d 2 x 

and 

J A R {\x\)d 2 x = J A{x)d 2 x, 
we get the desired property. If L is infinite, we can use a limiting procedure. Finally, we get 



J_ < 16tt 2 



d 2 xV R (\x\)in \x\ 



+ 32tt 2 J d 2 xV4x) J d 2 yV4y)(in(2+\y\)f 

A O 

J d 2 xV.{x) 



(A. 14) 



From (A. 2) and (A. 3) we see that J_ is bounded. This concludes the proof. 



One question is: can we go in the opposite direction? Assume that we know that (A.l) 
holds. There exists certainly a region \x — x a \ < R where Inf V~ = m > 0. If such a region did 
not exist, V~ would be zero almost everywhere! So 

/ > itR 2 m [ V~(x) £n[2 + \x\] 2 d 2 x (A. 15) 

J\x\>\x \+R+2 

Now we choose y > 4 + \x \ + R + 2 such that 

^\y-yo\<R' V -(v) =m' >0, 

then 

/ > irR' 2 m' ( V-(x)(£n(A - R')) 2 d 2 x. (A. 16) 

J\x\<\x \+R+2 

This proves that the convergence of (A.l) implies the convergence of (A. 2). 

It is not possible to deduce (A. 3) from (A.l) because (A. 3) involves Vr and (A.l) does not. 
However, in practice the conditions are very similar. Nevertheless, the following example shows 
that (A. 3) is stronger than (A.l), even for a potential which does not need rearrangement: take 
the central potential 

V ^ l) = ~\x\ 2 \£nx\ 2 (£n\£n\x\\^ iOT W < 2? 

= for \x\ > y~ ( A - 17 ) 

For 7 < | (A.l) and (A. 3) are divergent, for | < 7 < 1 (A.l) is convergent and (A. 3) is 
divergent, for 7 > 1 (A.l) and (A. 3) are convergent. 
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